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Abstract 

Multisymplectic systems, partial differential equations, fluid dynam- 
ics, conservation laws, potential vorticity We construct multisymplectic 
formulations of fluid dynamics using the inverse of the Lagrangian path 
map. This inverse map - the "back-to-labels" map - gives the initial 
Lagrangian label of the fluid particle that currently occupies each Eule- 
rian position. Explicitly enforcing the condition that the fluid particles 
carry their labels with the flow in Hamilton's principle leads to our mul- 
tisymplectic formulation. We use the multisymplectic one-form to obtain 
conservation laws for energy, momentum and an infinite set of conserva- 
tion laws arising from the particle-relabelling symmetry and leading to 
Kelvin's circulation theorem. We discuss how multisymplectic numerical 
integrators naturally arise in this approach. 



1 Introduction 

Definition 1. A system of partial differential equations (PDEs) is said to be 
multisymplectic if it is of the form 

K ai \ j _ dH 

where each of the two-forms 

K a = -K^ j (z)dz i Adz 1 

is closed. Here z is an ordered set of dependent variables, total differentiation 
with respect to each independent variable q a is denoted by the subscript a after 
a comma, and the Einstein summation convention is used. 

The closed two-form k° is associated with the independent variable q a \ it 
is analogous to the symplectic two-form for a Hamiltonian ordinary differential 
equation. Hence there is a symplectic structure associated with each indepen- 
dent variable. In the first of a series of papers, Bridges (1997) pioneered the 
development of multisymplectic systems, showing that the rich geometric struc- 
ture that is endowed by the symplectic two-forms can be used to understand 
the interaction and stability of nonlinear waves. For many important PDEs, the 
multisymplectic formulation has revealed hidden features that are important in 
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stability analysis. In order to preserve at least some of these features in numeri- 
cal simulations, Bridges & Reich (2001) introduced multisymplectic integrators, 
which generalise the symplectic methods that have been widely used in numeri- 
cal Hamiltonian dynamics. Hydon (2005) showed that multisymplectic systems 
of PDEs may be derived from Hamilton's principle whenever the Lagrangian 
is affine in the first-order derivatives and contains no higher-order derivatives. 
This can usually be achieved by introducing auxiliary variables to eliminate the 
derivatives. 

The aim of this paper is to provide a unified approach to producing multi- 
symplectic formulations of fluid dynamics, based on the inverse map. Our ap- 
proach covers all fluid dynamical equations that are written in Euler-Poincare 
form (Holm et al., 1998), i.e. all equations which arise due to the advection of 
fluid material. First we use the inverse map to form a canonical Eulcr-Lagrange 
equation (following the Clebsch representation given in Holm & Kupershmidt, 
1983). Then the Lagrangian is made affine in the space and time derivatives by 
using constraints that introduce additional variables. Following Hydon (2005), 
we obtain a one-form quasi-conservation law which, when it is pulled back to 
the space and time coordinates, gives conservation laws for momentum and en- 
ergy. We also obtain a two-form conservation law that represents conservation 
of symplecticity; when this is pulled back to the spatial coordinates, it leads to 
a conservation law for vorticity. The multisymplectic version of Noether's The- 
orem yields an infinite space of conservation laws from the particle-relabelling 
symmetry for fluid dynamics; these conservation laws imply Kelvin's circulation 
theorem. The conserved momentum that is canonically conjugate to the back- 
to-labels map plays a key role in the derivation of the conservation laws. The 
corresponding velocity is the convective velocity, whose geometric properties are 
discussed in Holm et al. (1986). 

In this paper we show how the above constructions are made in general, illus- 
trating this with examples. We also discuss how multisymplectic integrators can 
be constructed using these methods. Sections[S]and[3]review the relations among 
multisymplectic structures, the Clebsch representation and the momentum map 
associated with particle relabelling. Section [4] shows how to construct a multi- 
symplectic formulation of the Euler-Poincare equation for the diffcomorphism 
group (EPDiff), and derives the corresponding conservation laws, including the 
infinite set of conservation laws that yield Kelvin's circulation theorem. Sec- 
tion [5] extends this formulation to the Euler-Poincare equation with advected 
quantities. This is illustrated by the incompressible Euler equation, showing 
how the circulation theorem arises in the multisymplectic formulation. Section 
|B] sketches numerical issues in the multisymplectic framework. Finally, Section 
[7] summarises and outlines directions for future research. 

2 Review of multisymplectic structures 

This section reviews the formulation of multisymplectic systems and their con- 
servation laws, following Hydon (2005). 

A system of partial differential equation (PDEs) is multisymplectic provided 
that it can be represented as a variational problem with a Lagrangian that is 
affine in the first derivatives of the dependent variables: 

L = Lf(z)z%-H{z). (1) 
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The Euler-Lagrange equations are then 



dH 

where the functions 



KfMX^^v ( 2 ) 



K?Az) = -^r-^Pr (3) 

are coefficients of the multisymplectic structure matrix. We define the (closed) 
symplectic two-forms 

K a = ±K*{z)dz i Adzt, (4) 

and obtain the structural conservation law (Bridges, 1997). 

k% = 0. (5) 

Hydon showed that the Poincare Lemma leads to a one-form quasi-conservation 
law 

(LJdz*\ a = d(L°zi a - H{z)) = dL, (6) 

whose exterior derivative is 

Every one-parameter Lie group of point symmetries of the multisymplectic 
system ^ is generated by a differential operator of the form 

X = Q\q,z)^- + (Q\q,z)), a ^-. (7) 

Noether's Theorem implies that if X generates variational symmetries, that is, 
if 

XL = B% (8) 

for some functions B a , then the interior product of X with the one- form quasi- 
conservation law yields the conservation law 

(LjQ j - B a ), a = 0. (9) 

This is the multisymplectic form of Noether's theorem. 

Every multisymplectic system is invariant under translations in the inde- 
pendent variables q. For each of these symmetries, Noether's theorem yields a 
conservation law 

(Lfz^ ~ L5%), a = 0. 

Such conservation laws can equally well be obtained by pulling back the quasi- 
conservation law ([H]) to the base space of independent variables. Commonly, the 
independent variables are spatial position x and time t. Pulling back ([6|) to these 
base coordinates yields the energy conservation law from the dt component, and 
the momentum conservation law from the remaining components. We shall see 
the form of these conservation laws for fluid dynamics in later sections. 
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3 The inverse map and Clebsch representation 



3.1 Lagrangian fluid dynamics and the inverse map 

Lagrangian fluid dynamics provides evolution equations for particles moving 
with a fluid flow. This is typically done by writing down a flow map $ from 
some reference configuration to the fluid domain f2 at each instance in time. As 
the fluid particles cannot cavitate, superimpose or jump, this map must be a 
diffcomorphism. 

For an n-dimensional fluid flow, the flow map $ : R™ x R i— > R™ given 
by x = <&(/,£) specifies the spatial position at time t of the fluid particle that 
has label I = <I>(a;,0). The inverse map $~ 1 gives the label of the particle 
that occupies position x at time t as the function I = <E> -1 (a:,£). The Eulerian 
velocity field u(x, t) gives the velocity of the fluid particle that occupies position 
x at time t as follows: 

x(l, t) = u(x(l, t), t). 
Each label component lk(x,t) satisfies the advection law 

h,t + Uih,i = 0. (10) 

Here ^ and j2 ; denote differentiation with respect to t and x; respectively. We 
use Cartesian coordinates and the Euclidean inner produccl, so we shall not 
generally distinguish between 'up' and 'down' indices; summation from 1 to n 
is implied whenever an index is repeated. 

3.2 Clebsch representation using the inverse map 

A canonical variational principle for fluid dynamics may be formulated by follow- 
ing the standard Clebsch procedure using the inverse map (Seliger & Whitham 
(1968), Holm & Kupershmidt, 1983). The Clebsch procedure begins with a 
functional £[u] of the Eulerian fluid velocity «, which is known as the reduced 
Lagrangian in the context of Euler-Poincare reduction (Holm et al, 1998). One 
then enforces stationarity of the action S = J £[u]dt under the constraint that 
equation (|10p is satisfied by using a vector of n Lagrange multipliers, which 
is denoted as re. These Lagrange multipliers are the conjugate momenta to I 
in the course of the Legendre transformation to the Hamiltonian formulation. 
One may choose £[u] to be solely the kinetic energy, which depends only on 
u. More generally, t will also depend on thermodynamic Eulerian variables 
such as density, whose evolution may also be accommodated by introducing 
constraints. These constraints are often called the "Lin constraints" (Serrin, 
1959). This idea was also used in reformulating London's variational principle 
for superfluids (Lin, 1963). 

Definition 2 (Clebsch variational principle using inverse map). The Clebsch 
variational principle using the inverse map is 

S [ £[u]+ [ n-{l t +u-\'l)dV{x)dt = 0, 

1 This is only done for clarity and the equations are easily extended to the case when the 
domain Q is a curved manifold. 
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where iv{x,t) are Lagrange multipliers which enforce the constraint that particle 
labels l(x,t) are advected by the flow. 

Taking the indicated variations leads to the following equations: 

Su: + (V/) T ■ 7T - 0, (11) 

ou 

5tv : l t + (u- V)Z = 0, 
51: Tv t + V • (mr) = 0, 

where 

((VZ) T • tt) . := n k h,i, (V • (mr)) fe := (ujir^j, (12) 
and the variational derivative S£/Su is defined by 

f Sf 

t[u + eu'} = e\u] +e — u' 6V(x) + 0(e 2 ) . 
Jn° u 

Remark 3 (Clebsch representation). In the language of fluid mechanics, the 
expression ill]) for the spatial momentum m = 5£/Su in terms of canonically 
conjugate variables (Z,7r) is an example of a "Clebsch representation," which 
expresses the solution of the EPDiff equations (see below) in terms of canonical 
variables that evolve by standard canonical Hamilton equations. This has been 
known in the case of fluid mechanics for more than 100 years. For modern 
discussions of the Clebsch representation for ideal fluids, see, for example, Holm 
& Kupershmidt (1983) and Marsden & Weinstein (1983). In the language of 
geometric mechanics, the Clebsch representation is a momentum map. 

3.3 Particle relabelling 

As the physics of fluids should be independent of the labelling of particles, 
one may relabel the particles (by a diffeomorphism of the flow domain) with- 
out changing the dynamics. This is called the particle relabelling symmetry] 
Noether's theorem applied to this symmetry leads to the Kelvin circulation 
theorem. See Holm et al. (1998) for a modern description. 

3.4 Clebsch momentum map 

Definition 4. A momentum map is a map J : T*Q — > q* from the cotangent 
bundle T*Q of the configuration manifold Q to the dual g* of the Lie algebra q 
of a Lie group G that acts on Q. The momentum map is defined by the formula, 

%)-« = (",.6j(?)>, (13) 

where v q e T*Q and £ £ 0. In this formula £q is the infinitesimal generator of 
the action of G on Q associated with the Lie algebra element £, and {v q ,t;Q(q)} 
is the natural pairing of an element ofT*Q with an element ofT q Q. 

Theorem 5. The Clebsch relation ill]) defines a momentum map for the right 
action Diff (17) of the diffeomorphisms of the domain f2 on the back-to -labels map 

iE 

2 As we discuss later, this right action contrasts with fluid particle relabelling, which arises 
by the left action of the diffeomorphisms on the inverse map. 
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Proof. The spatial momentum in equation (jlip may be rewritten as a map 
Jo : T*fl i > 3C(fl) from the cotangent bundle of f2 to the dual X*(tt) of the 
vector fields X(f2) given by 

Jq : m • dec = - ((VZ) T • tt^J ■ dx = - tt ■ dl =: - 7r fe d^ . (14) 

That is, Jo maps the space of labels and their conjugate momenta (I, it) 6 T*fi 
to the space of one-form densities m £ X*(il) on il. The map (fT4"]) may be 
associated with the right action l-rj of smooth invertible maps (diffeomorphisms) 
r\ of the back-to-labels maps I by composition of functions, as follows, 

Diff(0) : l-r) = lor]. (15) 

The infinitesimal generator of this right action is obtained from its definition, 
as 

Xo(«):=-r (lo V (s))=TloX, (16) 



ds 



s=0 



in which the vector field X G 3£(f2) is tangent to the curve of diffeomorphisms rj s 
at the identity s = 0. Thus, pairing the map Jo with the vector field X G 3£(f2) 
yields 



(Jo(i,7T),X) =-(7T-dZ,X) 

= / 7r /c ^. J X i (a;)dV r (£c) 
Js 



= - ((l,ir),Tl>X) 

= -((!,*), x n (i)) , 

where ( • , • ) : Tj*f2 x XjSl M is the L 2 pairing of an element of Tj*f2 (a one- 
form density) with an element of TiQ (a vector field). 

Consequently, the Clebsch map (fTTj) satisfies the defining relation (fT5|) to be a 
momentum map, 

J(/,7r) = -7r-dZ, (17) 
with the L 2 pairing of the one- form density — 7r ■ dl with the vector field X. □ 

Remark 6. Being the cotangent lift of the action o/Diff(f2), the momentum 
map Jo in zs equivariant and Poisson. That is, substituting the canonical 
Poisson bracket into relation (|17[) yields the Lie-Poisson bracket on the space of 
m's. See, for example, Holm & Kupershmidt (1983) and Marsden & Weinstein 
(1983) for more explanation, discussion and applications. The momentum map 
property of the Clebsch representation guarantees that the canonically conjugate 
variables (Z, 7r) may be eliminated in favour of the spatial momentum m. Before 
its momentum map property was understood, the use of the Clebsch represen- 
tation to eliminate the canonical variables in favour of Eulerian fluid variables 
was a tantalising mystery (Seliger & Whitham, 1968). 

Note that the right action o/Diff(f2) on the inverse map is not a symmetry. 
In fact, as we shall see, the right action o/Diff(f2) on the inverse map generates 
the fluid motion itself. 
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3.5 Elimination theorem 

Eliminating the canonically conjugate variables (Z,7r) produces an equation of 
motion for m = 8£/8u, which is constructed in the proof of the following theo- 
rem: 

Theorem 7 (Elimination theorem). The labels I and their conjugate momenta 
7r may be eliminated from the equations arising from the variational principle 
JJ|) to obtain the weak form of the following equation of motion for S£/Su: 

d 81 81 

where 

ad* m — V • (um) + (Vit) T • m 

is defined by 

(ad^ m, w) = —(m, ad u w) = (m, (u ■ V)w — (w ■ V)it), 

and (•, •) is the L 2 inner-product. This equation of motion for 8£/Su is the 
Euler-P oincare equation for the diffeomorphism group (EPDiff) (Holm et al., 
1998). 

Proof. Take the time-derivative of the inner product of 5£/Su with a time- 
independent vector field w: 

UlA - 1/-™* -.-)-!( -,(.-vA 



V • (utt), (w ■ V)Z ) + ( 7T, (to ■ V)(u ■ V)Z ), 



= ^tt, -(it • V)(i£> • V)l + (w-V)(u- V)IJ = ^tt, - (ad„ w ■ V)Z 

= / - (VZ) T ■7r,ad„io\ = ^-^-,ad M u;^ = - ^ad* J^-,™^, 
which is the (weak form of the) EPDiff equation. □ 



3.6 Example: EPDiff^ 1 ) 

To give a concrete example, consider EPDiff with £[u] being the TJ^-norm for 
u. This is the n-dimensional Camassa-Holm (CH) equation (Camassa & Holm, 
1993; Holm et al., 1998; Holm & Marsden, 2004), which has applications in 
computational anatomy (Holm et al, 2004; Miller et al., 2002). This system 
has the reduced Lagrangian 

£[u]= f I (|w| 2 + A 2 lVw| 2 ) dV^Cas) = / I (u lUl + X 2 u^ UlJ ) dV(x) = l\\u\\ 2 Hl . 
Jn z Jn 1 1 A 

The EPDiff equation amounts to 

^ + (u ■ V)m + (Vm) t • m + mV • u = 0, m = (1 - A 2 V 2 )it . 
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When n = 1, these reduce to the Camassa-Holm (CH) equation, 



m t + um x + 2mu x = 0, 



TO = u — A It, 



3.7 Advected quantities 

To construct more general fluid equations we shall include advected quantities 
a whose flow-rules are defined by 



a t + C u a = 0, 



(18) 



where C u is the Lie derivative. Such advected variables typically arise in the 
potential energy or the thermodynamic internal energy of an ideal fluid. For 
example, advected scalars s (as in salinity) satisfy 



d_ 

dt 



s + C u s = 0, i.e. s f + (u ■ V)s = 0, 



and advected densities pdV satisfy. 
8 



dt 



( P dV)+C u (pdV) =0, i.e. Pt + y-(pu) = 0. 



(19) 



(20) 



A more extensive list of different types of advected quantity is given in Holm et 
al. (1998). 

We write the reduced Lagrangian £ as a functional of the Eulerian fluid 
variables u and a, and add further constraints to the action S to account for 
their advection relations, 

S = I i[u, a] dt + I dt I 7T • (l t + (u ■ V)Z) + <t>(a t + C u a) dV(x). (21) 
J J Jn 

The Euler-Lagrange equations, which follow from the stationarity condition 
SS = 0, are 



Su 

5tz 
SI 

5<t> 

Sa 



51 

— + (Vl) T --K + (j}oa = 0, 
ou 

l t + (u-V)l = 0, 

— 7T t — V • (UTT) = 0, 

at + C u a = 0, 

—<Pt - c u <f> + — = o, 

oa 



(22) 
(23) 
(24) 



where the diamond operator (o) is defined as the dual of the Lie derivative 
operation C u with respect to the L 2 pairing. Explicitly, under integration by 
parts, 

/ (<jx>a) -udV(x) = - [ ((f>£ u a)dV(x). (25) 



Remark 8. The map to the spatial momentum in equation 

81 



Su 



=: m = — 7r aVF 



cf> o a , 



(26) 
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is again a momentum map, this time for the semidirect-product action of the 
diffeomorphisms on Q x V* . Again the momentum map property allows the 
canonical variables to be eliminated in favour of the Eulerian quantities. As a 
result, eliminating the variables I, it and 4> leads to the Euler-Poincare equation 
with advected quantities a. 

Theorem 9 (Elimination theorem with advected quantities). The labels I, their 
conjugate momenta ir and the conjugate momentum (</>) to the advected quan- 
tities (a) may be eliminated from equations 122\24\) to obtain the weak form of 
the Euler-Poincare equation with advected quantities: 

d si si 5i r 

— - had„— = ao— , a t + £ u a = Q. 

at ou ou oa 

Proof. Take the time-derivative of the inner product of 6£/6u with a function 
of w: 

jtitA - U- mT ^-^ w ) = U- w ' {w ^ )i ) + U^ Cwa \ 



= • (mr), (w ■ V)l\ + ^tt, (w ■ V)(tt • V)l\ + p«a,w 

+/ — ^ - C w a^j + (^(j), -C w C u a^j , 
= lir,-(u • V)(u> • V)l + (w ■ V)(u ■ V)z\ + / y oa,w\ 



= ^tt, - (ad u to • V)Z^ + /^,£ad ul0 a\ + |^oa,«)|, 

= / - (VZ) T • 7r - (jx>a, ad„M)^ + oa.ic|, 

/ Si \ / Si \ I ^ 5£ 51 \ 

□ 

Remark 10. These Euler-Poincare equations with advected quantities cover 
all conservative fluid equations which describe the advection of material. For a 
large collection of examples, see Holm ct al. (1998). 

3.8 Example: Incompressible Euler equations 

As an example, consider the reduced Lagrangian for the incompressible Euler 
equations 



l[u,p,p}= J^^L+p(l-p)dV(x). 
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Here p(x, t) is the ratio of the local fluid density to the average density over f2; 
this is governed by the continuity equation (j2T)|) . The pressure p is a Lagrange 
multiplier that fixes the incompressibility constraint p = 1. The variational 
derivatives in this case are 

61 _ Si _ \u\ 2 81 _ . 

Su ^ ' Sp 2 ^' Sp ^ 

Consequently, the Euler-Poincare equations become 

(pu) t + (u ■ V)(pu) + pu(V ■ u) + p{Vu) T ■ u = pV 

p t + V • (pu) = 0, 
P = 1, 

and rearrangement gives the Euler fluid equations, 

u t + (u ■ V)u = -Vp, V • u = 0. 



4 Inverse map multisymplectic formulation for 
EPDiff(// 1 ) 

As we now have a canonical variational principle for fluid dynamics via the 
inverse map, one may obtain its multisymplectic formulation by extending the 
phase space so that the Lagrangian is affine in the space and time derivatives. In 
this section we show how to do this for EPDiff(i/ 1 ) as discussed in the previous 
section. 

4.1 Affine Lagrangian for EPDiff(if 1 ) 

After introducing the inverse map constraint, the Lagrangian becomes 

1 A 2 

L = -UiUi + —Uiju^j + n k {l k<t + Ujlk,j) ■ 

Any high-order derivatives and nonlinear functions of first-order derivatives 
must now be removed from the Lagrangian to make it affine. We introduce 
a tensor variable 

Wij = u itj ; 

this relationship may be enforced by using Lagrange multipliers. However, it 
turns out that the multipliers can be eliminated and the Lagrangian becomes 

1 A 2 

L = -UtUi - —WijWij + A J lt ,,</,., + TT k (l k>t + Ujl k ,j) , (27) 
which is now affine in the space and time derivatives of u, W, I and tt. 
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4.2 Multisymplectic structure 

The Euler-Lagrange equations for the affine Lagrangian (T2T|) are 



Sui 

5n k 
SWu 



Ui - A"' II ', ,. , +7Tfcife,i = 0, 
-TTfc.i - (TTfcUj)j = 0. 

-A 2 W l? - + A 2 u,-, i = 0. 



These equations possess the following multisymplectic structure as in equation 



/ TTfcft 




-a 2 ^ 






-TTfe^ -9 t 











dt + Ujdj 












V A 2 a, o 





o J 






= d/dt, di = d/dxi, and 










/l A 2 








X 2 


i? = - ( -nm - —Wij 


= 


-Gh. 


2 


4"' 



4.3 One-form quasi-conservation law 

For our multisymplectic formulation of EPDiff(i7 1 ), the independent variables 
are 



q J =x.j, j = l, 



t- 



and the dependent variables are 



n-\-k 



z 2n+k = n k , 



y (i+2)n+j = w 



(28) 



where i,j and k range from 1 to n. Comparing (|27p with ([T]) gives the following 
non-zero components L": 



,,72. 



L\ = X 2 Wij , L J n+k = TT kUj , LIU = n k , i,j,k=l, 
Therefore the one-form quasi-conservation law amounts to 

(7r fc di fc ) it + (X 2 Wijdui + n k Ujdlk) d = dL. (29) 
The exterior derivative of this expression yields the structural conservation law 

(dn k A dl k ) t + (\ 2 dWij A du, + Ujdir k A dl k + n k duj A dl k ) = 0. (30) 

4.4 Conservation of energy 

For EPDiff(i? 1 ), the di-component of the pullback of the one-form conservation 
law flU gives 

(Tfkh.t ~ L) t + (X 2 W tJ Ui tt + TT k Ujl k ,t) t 
In terms of u and its derivatives, this amounts to 
1 A 2 
2 



0. 



™ " 2 UiUi 



-UijUijj + (\ 2 UijUij + UiUjirii) . = 0, 



where 



rrii = Ui — X Ui 



kk ■ 



This is the energy conservation law for EPDiff(i7 1 ). 
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4.5 Conservation of momentum 

Similarly, the conservation law that is associated with translations in the Xi 
direction is 

{■7Tklk,i),t + (A 2 WfejUfe,i + TTkUjh,i - SijL) . = 0, 

which amounts to the momentum conservation law 

4.6 Conservation of vorticity 

Next, consider the coefficient of each dx r A da; s in the pull-back of the structural 
(two-form) conservation law (j5U|) . This is 

(ftk,rlk,s — 7Tfe,s^fe,r) t + (A 2 (Wy ir « ijS — Wij, s Ui >r ) + Uj {lTk,rlk,s - TTk, s lk,r) + ^k {Uj, r lk 

which amounts to 

(m r , s - m s , r ) t + (A 2 (Ui, s Ui t j r - Ui^Uijs) + (ujUlr)^ - (uj-m s ), r ) . = 0. 

One can regard this as a vorticity conservation law for EPDiff(iJ 1 ); it is a 
differential consequence of the momentum conservation law. 

4.7 Particle relabelling symmetry 

As we discussed in Section[3j fluid equations in general, and EPDiff in particular, 
are invariant under relabelling of particles. In the context of the inverse map 
variables, relabelling is accomplished by the action of the diffeomorphism group 
Diff (fi) defined by 

i^7]oi = n(i), ?/ g Diff(n). 

The corresponding infinitesimal action of the vector fields 3£(Q) is then 

Z^£oi = £(Z), ^I(O), 
and the cotangent lift of this action is 

(7r,0^(-(V£(0) T .7r,£(0). 

To obtain the symmetry generator 0, we extend the above action to first 
derivatives as follows: 



dlk aht ol 



k . i 



QUI) d ( dOS)\ d ( d$ k {l)\ d 

-TTk—^: T{ n k—^ T{ TTfc- 



dlj dnj \ dlj J , diTjj V ®h ) i 
The relabelling symmetries are variational, because 

XL = nk(£k(l)),t + KkUi(£k(l)),i - 7r ^ ~^~" + u ih,i) = °- 
Noether's theorem then gives the conservation law 

(*-k&(0),t + ^kuMi)),, = o. 

A conservation law exists for each element £ of X(f2), so particle relabelling 
generates an infinite space of conservation laws. 
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4.8 Circulation theorem 

To see how the particle relabelling conservation laws relate to conservation of 
circulation, note that if p is any density that satisfies 

P,t + {pu>j),j = 0, 

then 

^MT) + Uj ( 7r * e * (0 ) =o. 

If we pick a loop C(t) which is advected with the flow, then 

dt Jc(t) P 

For a vector field £ which is tangent to the loop at time 0, and satisfies |£| = 1 
on the loop, then 

£dx = (VI) ■ da; 

for all times t, and one finds 

U JLM.d-O. (3!) 
dt Jc(t) P 

The momentum formula (|11[) gives 

d r (i-a 2 v 2 )m , d r m , 

' • dec = — <z> — • da; = 0, 



/c(t) P dt J c[t) p 

which is the circulation theorem for EPDiff. 

5 Inverse map multisymplectic formulation for 
Euler-Poincare equation with advected quan- 
tities 

To extend this method to more general equations with advected quantities is 
very simple: take the Lagrangian obtained from equation (|3|) and add variables 
to represent higher-order derivatives. For the sake of brevity we shall com- 
pute one example, the incompressible Euler equations, and briefly discuss the 
implications for the circulation theorem. 

5.1 Multisymplectic form of incompressible Euler equa- 
tions 

We start with the reduced Lagrangian 

t[u,p,p]= I \puiUi +p(l - p)dV(x), 
Jn z 

where p is the pressure and p is the relative density, and add dynamical con- 
straints to form the Lagrangian: 

L = ^puiUi + p(l - p) + n k (l kit + Uilk,i) + 4> (p,t + (pui),i) ■ 
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This Lagrangian is already affine in the first-order derivatives, so the Euler- 
Lagrange equations are automatically multisymplectic in these variables: 
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where the quantity 

H = - (^puiu t +p(l - p) 
is negative of the Hamiltonian density. 



5.2 Circulation theorem for advected quantities 

The conservation law for particle-relabelling follows exactly as in Section [H and 
we obtain equation ()31|) as before. The difference is that now the momentum 
formula (momentum map) is 

m — — = -n k V/fe - <p o a 
du 

and so one obtains 

d f m , f ld£ 

— <p — ■ ax = <p — — o a ■ ax. 

dt J C (t) P Jc(t) P oa 

For the incompressible Euler equations, a is the relative density p, so 

d£ ^d£ 
— oa = pV— , 
da dp 

which leads to the circulation theorem 

d I m , I r-,M , 

— <p — • ax = <t> V— • da; = 0. 

dtJ C {t) P Jc(t) <JP 



6 A note on multisymplectic integrators 

In this section we discuss briefly how to produce multisymplectic numerical 
integrators, using the inverse map formulation given in this paper. We note 
in particular that the multisymplectic method will satisfy a discrete form of 
the particle-relabelling symmetry and hence we will obtain a method that has 
discrete conservation laws for — 7r • VZ. 



6.1 Variational integrators 

A multisymplectic integrator for a PDE is a numerical method which preserves 
a discrete conservation law for the two-form k given in equation (Bridges 
& Reich, 2001). As described in (Hydon, 2005), a discrete variational principle 
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with a Lagrangian that is affinc in first-order differences automatically leads 
to a set of difference equations which are multisymplectic. This now makes it 
very simple to construct multisymplectic integrators for fluid dynamics using 
the inverse map formulation: one simply replaces the spatial and time integrals 
in the action with numerical quadratures, replaces the first-order derivatives 
by differences, and takes variations following the standard variational integra- 
tor approach (Lew et ai, 2003). Whilst the method will preserve the discrete 
conservation law for the two-form k, the one-form quasi-conservation law will 
not be preserved in general, and hence the other conservation laws will not be 
exactly preserved. 

6.2 Discrete relabelling symmetry 

As 7T and I are still continuous in the discretised equations, the multisymplectic 
integrator will have a discrete particle-relabelling symmetry analogous to the 
one given in Section IU with the only difference being the discretisation of the 
cotangent lift. Following the variational integrator programme described in Lew 
et at (2003), the discrete form of Noether's theorem will give rise to discrete 
conservation laws for the multisymplectic method. 

6.3 Remapping labels 

If this approach is to be applied to numerical solutions with intense vorticity then 
one needs to address the problem that eventually the numerical discretisation of 
the labels I will become very poor due to tangling, and hence the approximation 
to the momentum 

m = -(Vl) T ir = -ir k Vl kl (32) 

will degrade with time. One possible approach would be to apply discrete 
particle-relabelling, mapping the labels back to the Eulerian grid in such a 
way that the momentum (|32[) stays fixed. This transformation is exactly the 
relabelling given in Sectional Numerically, one could construct a transformation 
(using a generating function for example) which satisfies 

I h-> X + 0(Ax p ,At p ), tt(VJ) i — * 7r(VZ) + 0{Ax p , At p ), 

where p is the order of the method. For instance, one might use a variational 
discretisation of the relabelling transformation, which is generated by a sym- 
plectic vector field whose Hamiltonian is tt ■ £(Z) = 7Tfe£fc(i). In this way, one 
may still retain some of the conservative properties of the method. 

7 Summary and Outlook 
7.1 Summary 

This paper describes a multisymplectic formulation of Eulcr-Poincare equations 
(which arc, in essence, fluid dynamical equations with a particle-relabelling 
symmetry). We have used the inverse map to obtain a canonical variational 
principle, following Holm and Kupershmidt (1983). As noted in Hydon (2005), 
a multisymplectic formulation can be obtained by choosing variables such that 
the Lagrangian at most linear in the first-order derivatives, and contains no 
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higher-order derivatives. We have shown how to construct the multisymplectic 
formulation for the Euler-Poincare equations for diffcomorphisms, using the 
example of the EPDiff(i/ 1 ) equations, and how to extend the method to the 
Euler-Poincare equations with advected quantities. These equations encompass 
many fluid systems, including incompressible Eulcr, shallow-water, Euler-alpha, 
Green-Naghdi, perfect complex fluids, inviscid magnctohydrodynamics, etc. 

The techniques of Hydon (2005) have led to conservation laws for these 
systems, including the usual multisymplectic conservation laws for energy and 
momentum plus an infinite set of conservation laws which arise from the particle- 
relabelling symmetry of fluid dynamics. We have highlighted the connection 
between these latter conservation laws and Kelvin's circulation theorem, and 
showed that multisymplectic integrators based on this formulation will have 
discrete conservation laws associated with this symmetry. 

7.2 Outlook 

In the last section of this paper we have discussed the possibility of developing 
multisymplectic integrators for fluids using this framework. It is undoubtedly 
simple to construct such integrators, but the issue of accuracy with time arises 
whenever the flow is strongly mixing and numerical errors make the label field I 
very noisy. A discretisation of the relabelling map discussed in this paper could 
provide a way to prevent this problem whilst retaining some of the geometric 
properties of the method. These ideas may aid the future development of in- 
tegrators that have conservation laws for vorticity and circulation, which are 
desirable for numerical weather prediction and other applications. 

In a different direction, we believe that multisymplectic integrators would be 
especially apt for applications of EPDiff to template-matching in computational 
anatomy (Holm et al, 2004). The matching problem is an initial- final value 
problem. In such problems, space and time may be treated on an equal footing, 
just as in the multisymplectic formulation. 
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